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ABSTRACT 



PSR J1910-5959A is a binary pulsar with a helium white dwarf compan- 
ion located about 6 arcmin from the center of the globular cluster NGC6752. 
Based on 12 years of observations at the Parkes radio telescope, the relativistic 
Shapiro delay has been detected in this system. We obtain a companion mass 
Mq = 0.180 ± O.O18M (la) implying that the pulsar mass lies in the range 
1.1 M Q < Mp < 1.5 Mq. We compare our results with previous optical de- 
terminations of the companion mass, and examine prospects for using this new 
measurement for calibrating the mass— radius relation for helium white dwarfs 
and for investigating their evolution in a pulsar binary system. Finally we ex- 
amine the set of binary systems hosting a millisecond pulsar and a low-mass 
helium white dwarf for which the mass of both stars has been measured. We 
confirm that the correlation bet ween the companion mass and the orbital period 



predicted by iTauris fc Savonijel reproduces the observed values but find that the 
predicted Mp — P B correlation over-estimates the neutron-star mass by about 
O.5M in the orbital period range covered by the observations. Moreover, a few 
systems do not obey the observed M P — P B correlation. We discuss these results 
in the framework of the mechanisms that inhibit the accretion of matter by a 
neutron star during its evolution in a low-mass X-ray binary. 



Subject headings: pulsars: individual (PSR J1910— 5959A) — globular clusters: 
individual (NGC 6752) 



-3 - 



INTRODUCTION 



PSR J1910-5959A (henceforth PSR A) is a millisecond pu lsar (MSP) discovered in the 



Parkes Globular Cluster Pulsar Survey (jD'Amico et al. 



2001 



hereafter paper I). It has a 



spin period of 3.27ms, and it orbits around a low-mass companion with Mq > 0.2 M Q , 
assuming a pulsar mass Mp = 1.4 M , with an orbital period Pb = 0.84 days. Because of 
its unusual position with respect to the core of the globular cluster NGC 6752 (its angular 
distance from the cluster center is 8± = 6'A, corresponding to 3.3 half- mass radii; paper I), 



it has been questioned whether this object is a member of the cluster ( IBassa et al. 
In fact, of all cluster-; 
of the hosting cluster 



2006|). 



In fact, of all cluster-associated MSPs, PSRA has the largest angular distance from the core 



The bi nary companion of P SRA has been recogni sed as a helium w hite dwarf star 



(HeWD) by 



Ferraro et al. 



(120031 ) and independently by 



Bassa et al. 



( 120031 ). via observations 



taken with the ESO Very Large Telescope (VLT) and the Hubble Space Telescope (HST). 
The values obtained by these authors for the mass of the binary companion of PSRA 
are in agreement (M co ~ 0.17 — 0.20 M Q ). Both authors also agree in concluding that 
the photometric properties of this star are consistent with the hypothesis that PSRA is 
associated with the globular cluster NGC 6752. 

The radia l veloc ity c urve of the compan ion of PSRA has been investigated by 



Cocozza et al 



(|2006h and 



Bassa et al. 



( 120061 ) . through phase- resolved spectroscopy in 



the optical band with HST observations. The resulting values for the mass ratio of the 
binary system, for the mass of the optical companion and for the orbital inclination were 



also compatible with each other, although 



Bassa et al 



(120061 ) obtained more stringent 



Tntriguingly PSR J1910-5959C in the same cluster ranks second in this list, being located 
at 9± = 2'.7, corresponding to 1.4 half-mass radii, from the core of the cluster (paper I). 
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constraints (see §3.21 for details). iBassa et al.l (120061 ) carefully discussed the association of 
this binary system with NGC 6752, and indicated a preference for non-association, though 
their arguments are not conclusive because of the large uncertainties on the parameters 
that would allow one to discriminate between the two scenarios (see §3.21 for details). 



Corongiu et al. 



(120061) (hereafter 



A refined ephemeris for PSRA was presented by 
paper II), where a low eccentricity model was applied to describe the pulsar's orbit. The 
measured eccentricity was e ~ 3 x 10~ 6 at a longitude of periastron u ~ 90°. Paper II 
discussed the possibility that such eccentricity was indeed a misinterpreted Shapiro delay. 
Nevertheless, paper II maintained a conservative approach to this problem and deferred the 
revision of this issue until a significantly longer time span of timing data became available. 

Correlations betwe en the orbital peri o d, th e pulsar mass and the companion mass 



have been obtained by 



Tauris fc Savonijd (119991 . hereafter TS99), from their numerical 



calculations on the non-conservative evolution of close binary systems with low-mass donor 
stars (Me < 2M ), a 1.3 M & accreting neutron star and an orbital period Pp > 2 days 
at the beginning of their low-mass X-ray binary (LMXB) phase. These systems have 



diverging or 
(LMBMSP. 



jital separation and are the progenito rs of low-mass binary millisecond pulsars 



Bhattacharya fc van den Heuvellll99ll ). TS99 obtained a positive correlation 



between the orbital period and the companion mass and a negative correlation between the 
or bital period and the pulsar mass . The pr e dicted Mq — Pb correlation has been observed 



by 



van Kerkwijk et al. 



(J2005J) and 



Lorimerl (120081 ). TS99 remarked on the discrepancy 



between the predicted Mp — Pb correlation and the available data. However in most cases 
the determinations of the Mq and Mp were based on models for the companion atmosphere 
or on statistical hypothesis for the orbital inclination, often leading to large uncertainties. 
The observation of a Shapiro delay allows one to obtain a direct measurement of Mc and 
Mp and gives a strong test for the aforementioned correlations. 
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In this paper we present timing data for PSRA collected over 12 years at the Parkes 
radio telescope. This data span is twice than that used in paper II and allowed us to get 
better constraints on the rotational, positional and orbital parameters for PSRA, and to 
detect the Shapiro delay in this system. The paper is organized as follows. In section §[2] 
we briefly describe the observing methods and present the timing solution based on the 
available data span, in §3] we describe the method for measuring the Shapiro delay, present 
our determination of the companion mass and discuss our results and their implications. In 
§4] we summarize the measurements for the masses of MSPs and their HeWD companions, 
focussing on those binaries for which Mq and Mp have been obtained from the detection 
of the Shapiro delay effect. We then compare these data to the predictions of TS99 and 
discuss the results. In $5] we briefly summarise our work. 



2. OBSERVATIONS AND GENERAL TIMING RESULTS 



Regular pulsar timing observations of PSRA in NGC 6752 have been carried out since 



September 2000 with the Parkes 64 m radio telescope at a central frequency of 



390 M Hz, 



1996|) 



using the central beam of the 20-cm multibeam receiver (jStaveley-Smith et al. 
or the H-OH receiver. The hardware system is the same as that used in the discovery 
observations (paper I). The effects of interstellar dispersion are minimised by using a 
filterbank having 512x0.5 MHz frequency channels for each polarisation. After detection, 
the signals from individual channels are added in polarisation pairs, integrated, 1-bit 
digitised every 80 /xs (125 /xs in earlier observations), and recorded to magnetic tape for 
off-line analysis. We synchronously folded the data at the pulsar period using the program 



DSPSR (Ivan Straten fc Bailes 



20111 ) with a subintegration time of 1 min and dedispersion in 



groups of 64 channels to give eight frequency sub-bands. We visually inspected each file to 
extract the maximum number of pulse times of arrivals (ToAs), depending on integration 
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time and interstellar scintillation which causes large fluctuations in the detected flux density. 
We considered a ToA reliable if the pulse profile had a signal-to-noise ratio greater than 
5 and the pulse was visible. This approach gave us ~ 1000 ToAs at different frequencies 
within the band of the receivers and consequently allowed us to fit for the dispersion 
measure and its first derivative. We calculated the ToAs by fitting a template profile to the 
observed mean pulse profiles and wc analysed them using the program TEMP02 (Hobbs et al. 



20061 ) with the DE405 solar system ephemeris and the TT(TAI) reference timescale. Tabled] 
summarises the best fit values and la uncertainties from TEMP02 (everywhere in the paper 
all reported measures are quoted with their la uncertainty, and all reported ranges are at 
the same confidence level). The timing residuals for the fit which have a \ 2 of 989.3 (with 
989 degrees of freedom) are shown in Fig.[TJ against MJD in the upper panel and against 
orbital phase in the lower panel. We determined the values and the uncertainties for the 
Shapiro delay parameters through a Bayesian analysis (see $3]). 

The new positional and rotational parameters at the reference epoch are all consistent 
with those reported in paper I and paper II at a confidence level of 2a, with the only 
exception being the proper motion, whose components in right ascension and declination 
are consistent at the 4a and 6a level respectively. The value for the first derivative of the 
dispersion measure (DM1) is of the same order of magnitude of similar determinations for 
other MSPs. The consistency of the orbital parameters will be discussed in §31 

3. THE SHAPIRO DELAY MEASUREMENT 
3.1. Method 



The orbital solutions presented and discussed in papers I and II were based on 
the ELL1 binary model, which is suitable for systems with small orbital eccentricities 
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(Wex 1998 unpublished work; see lLange et al.ll200ll ). In paper II an orbital eccentricity 
e = (3.4 ± 0.6) x 10~ 6 was reported, at a longitude of periastron u = 90° ± 10°. The 
values obtained for e and in particular for u were already commented as a possible 
misinterpretation of a Shapiro delay (paper II). By applying the same binary model of 
paper II to the longer data span presented in this work, we obtain e = (3.5 ± 0.5) x 10 -6 
at a longitude of periastron u = 95° ± 7°. Since the uncertainties on these two parameters 
did not decrease as expected for the case of a pure Keplerian ELL1 model we carefully 
reinvestigated the possible presence of a Shapiro Delay. 



We used the DD binary model (IDamour fc Deruelle 



1985 



Damour fc Deruelle 



19861 ) 



and included the two parameters that describe the Shapiro delay, namely the range r and 
the shape s. These two parameters are directly related to the companion mass Mq and the 
orbital inclination i through the simple relations r = T Q Mc (T & = GM Q /c 3 , where G 
is N ewton's universal gravitat ion constant and c is the speed of light) and s = sini (see, 



e.g. 



Lorimer fc Kramer 



20051 ). Fig. [2] displays our attempts to detect the first and third 



harm onic of the Shapiro delay in our timing residuals. We recall (e.g. 



Lorimer fc Kramer 



20051 ) that the first harmonic is obtained by fitting for all parameters, including the ones 
describing the Shapiro delay, then resetting Mc = 0.0 and sin 2 = 1.0 while keeping 
all other parameters to their best fit values, and finally by plotting the resulting timing 
residuals. The third harmonic is obtained by fitting for all parameters but the Shapiro 
delay ones, which remain set to the values Mc = 0.0 and sini = 1.0, then plotting the 
residuals for this best fit solution obtained in this way. While the first harmonic is still 
consistent with the response for an elliptic orbit, the third harmonic is a unique signature 
of the Shapiro delay. In the upper panels of Fig.[2] timing residuals show some evidence of 
the harmonic structure, while in lower panels, where residuals have been rebinned in 20 
orbital phase bins, the two harmonics are clearly recognisable. As a further comparison, in 
the lower panels we plot (solid lines) the theoretical curves for the Shapiro delay with which 
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our rebinned residuals are in satisfactory agreement. 



As an additional sanity check we have checked the consistency of the three Keplerian 
parameters that are common to the two different orbital models presented in paper II 
and in this work. The values for Pb are fully consistent to each other, and the same 
is true comparing T^sc m paper II and T here. We recall that in the case of zero 
eccentricity orbits, the periastron is defined as coincident to the ascending node. The 
only parameter for which we now find an inconsistency at a level of about 10a with 
respect to paper II is the orbit projected semi-mayor axis ap sin i. However because of 
the change in the adopted binary model (from ELL1 to DD) this parameter cannot 
be directly compared. The quantity to be compared between is ap sin z/(l + e). From 
paper II we obtain ap sin«/(l + e) = 1.2060425 ± 0.0000007, while in this work it results 
ap sini = 1.2060420 ± 0.0000008, as the orbital eccentricity has been now set to zero. The 
two values are then in full agreement. 

We obtained the values for the companion mass and the inclinat ion angle through a 



Bayesian analysis on our timing residuals. The method is described in 
We chose the a priori probability density functions as follows: 



Splaver et al. 



f|2002h . 



a) A flat distribution for cosi in the range < cosi < 1, as we assumed a randomly 
oriented orbit in space 



b) A flat distribution for Mq in the range O.OM < Mq < 0.5 M , chosen for 
consistency with the values of th is parameter resu l ting from the study of the radial 



velocity curve of the companion (IBassa et al. 



2006 



Cocozza et al. 



20061 ) 



Fig. [3] displays in panel a) the A% 2 map obtained after calculating the x 2 values in 
a 2048x2048 grid for the parameters of our interest. In panels b) and c) the posterior 
probability density functions for cos i and Mq respectively are reported. From this analysis 
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we derived a lower limit for the orbital inclination sini > 0.9994, i.e. % > 88°, and the 
value for the companion mass Mq = 0.180 ± 0.018 M@. 



3.2. Implications 



Under the hypothesis that PSRA is associated with NGC6752, the measurement 
reported in §3.11 is the second of this kind for a bi nary pulsar in a globular cluster after 



PSRJ1807-2459B in NGC6544 (ILvnch et al 



20121 ). PSRJ1807-2459B is a millisecond 



pulsar with a spin period Pq = 4.19 ms in a highly eccentric orbit e = 0.747. The nature of 
the companion is not yet determined, however its high mass Mq = 1.21 M & points towards 
a neutron star or a massive white dwarf companion. Thus PSRA represents the first case 
of a low-mass binary pulsar in a globular cluster for which the Shapiro delay has been 
detected. 

Our determination of the companion mass, Mq = 0.180 ± 0.018 M , is solely based 
on the relativistic Shapiro delay and so we did not assume any a priori hypothesis on 
the nature of the companion or the modelling of its structure. In this respect the full 
agreement with the values for Mq obtained using observations in the optical band provides 



independent support to the previously suggested interpretation of the na ture of the 



companion of PSRA, i.e., that it is a 



Bassa et al 



2006 



Ferraro et al. 



2003 



ow-mass helium-cor e white dwarf (IBassa et al 



Cocozza et al 



jmary 



2003 . 



20061 ). This also gives support to the 



physical modelling (i.e. the relationship between mass, radius, temperature and surface 
gravity of this category of white dwarfs) invoked for deriving the parameters of PSRA from 
the photometric and spectroscopic data only. 



Bassa et al. 



(120061 ) attempted to use a comparison between the observed parameters and 



the models mentioned above in order to decide whether or not this binary system belongs 



- 10 - 



to the globular cluster NGC6752. Unfortunately the uncertainties in the measurements 
of the companion parameters prevented the m from drawing; a firm conclusion. Since 



our uncertainty on Mq is similar to that of 



Bassaetal 



(120061 ) . our result does not help 



in addressing this point. Further support for the association of PSRA to NGC6752 is 
expected from the measurement of proper motion for the three isolated pulsars in the core 
of NGC 6752^ and comparison with the already accurately determined proper motion for 
PSRA (paper II). 

Assuming that PSRA is associated with NGC 6752, the timing determination of 
the companion mass opens the possibility of using the system for testing the various 
theoretical mass-radius relations proposed for helium white dwarfs. The large uncertainties 
in the determination of the masses and of the distances of the white dwarfs via optical 
observations make this difficult. This problem is particularly important for the low-mass 
white dwarfs, for which accurate determinations of the mass and the radius are rare. In 
this case we have independent measurements for the mass of the companion of PSRA 
and its distance, which is equal to the cluster's one under our hypothesis. From optical 



data 



BassaetaL 



(120061 ) derive a companion radius, namely Rq = 0.058 ± 0.004 R Q if 



PSRA belong s to NGC 6752. We not e tha t the mass-radi u s rela tions predicted by both 



the models of 



Rohrmann et al. 



(120021 ) and 



Serenelli et al. 



( 120021 ) nicely agree with the 



case of a white dwarf of R c = 0.058 ± 0.004 i? and M n = 0.180 ± 0.0 1 8 The two 



models differ in the progenitor's metallicity, namely the 
low-mass white dwarfs with solar metallicity while the 



Rohrmann et a 



Serenelli et al. 



(12002 



(I2002h describes 



model is 



suitable for metallicities comparable to the one of NGC 6752. The model by iDriebe et al. 



2 They certainly belong to NGC 6752 since at least two of them show the large \P\ at- 
tributed to the gravitational pull of the cluster. Their proper motion has not been measured 
yet, since they are weaker sources. 
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( 119981 ) computes the white dwarf radiu s only for masses Mwn > 0.180 M . It cannot be 



yet discarded (see for reference fig. 6 in 



Bassaetal 



20061 ) 



Given a typical ~8% uncertainty in the radius determination from optical observations, 
in order to discriminate between various models for the mass-radius relations, an uncertainty 
on Mq of 0.002 M Q is typically required. Unfortunately several decades are required to 
achieve this goal with a 64-m class telescope; such high precision could be obtained in about 
one decade by observing PSRA with the Square Kilometer Array. 

A determination of the mass Mp for PSRA can be obtained from our Shapiro delay 
measurement. From the lower limit in the orbital inclination and the pulsar mass function, 
we can obtain a firm range for Mp, 1.1 M & < Mp < 1.5 M , from which we deduce the 
total mass for the system lying in the range 1.26 M < M tot < 1.70 M & . If we instead use 



compa nion' s mass function 



(120061 ) and 



Cocozza et al. 



rom the studies on the companion's velocity curve by 



(120061 ). we can slightly improve our results. 



curve by 


Bassa et al. 


Bassa et al. 


( 


2006 


) 



report Mp/Mc = 7.36 ± 0.25, from which we now obtain M tot = 1.50 ± 0.14 M and 



Cocozza et al. 



(120061 ) report M P /M C = 7.49 ± 0.68, from 



M P = 1.32 ± 0.14 M , while 
which we now obtain M tot = 1.53 ± 0.18 M and M P = 1.35 ± 0.18 M . The weighted 
average of these two values yields Mp = 1.33 ± 0.11 M . 

It is worth mentioning that the value for the mass of PSRA is close to the lower edge of 
the mass range for binary millisecond pulsars located in the Galactic disk whose companion 
is a He WD star. Since in these systems the neutron star has experienced only one recycling 
phase, thus leaving the original neutron star mass almost unchanged, the value for the mass 
of PSRA suggests that this pulsar underwent only one recycling stage. If we combine this 



statement with the hypothesis that PSRA belongs to t 



re cluster, we can 



ejection of this system towards the cluster's outskirts ( IColpi et al. 



2002 



speculate that th e 



Colpi et al. 



2003) 
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most likely happened in an early stage of the cluster evolution, preventing this binary from 
experiencing further exchange interactions with the stars in the crowded environment of the 
cluster's core. This appears to be at odds with the predictions o f the models which suggest 



that the ejection occurred in the last ~ 1 Gyr (jColpi et al. 



20021). 



4. CORRELATIONS AMONG THE ORBITAL PERIOD, THE PULSAR 
AND THE COMPANION MASS FOR LOW MASS BINARY PULSARS 

TS99 presented detailed numerical calculations on the non- conservative evolution of 
close binary systems with a low-mass (Md < 2.0 M©) donor star, a 1.3 M & accreting 
neutron star and an orbital period Pb > 2 days at the b eginning of the LMXB ph ase. 



The evolution of LMXBs has already been modelled by 



Pylyser fc Savonijel (11988 



19891 ). who predict that binaries with starting periods Pb > Puf — 2 days increase 
their orbital separation, hence their orbital period, during the mass transfer (diverging 
systems), while the ones with Pg < Pbif decrease their orbital separation (converging 
systems). In this second case, a simple argument shows that the binary period decreases if 
\M tot \/M tot < \a\/a and increases if \M tot \/M tot > \a\/a, where M tot is the total mass of 
the binary and a is its orbital separation. 

TS99 focussed their attention on diverging LMXBs and they obtain a positive 
correlation between the orbital period and the final mass of the He WD companion, i.e., 
its mass increases with the observed orbital period, and a negative correlation between 
the orbital period and the mass of the pulsar. In order to compare the theoretical results 
of TS99 to the observed systems, we searched the literature for low-mass binary pulsars 
for which a reliable measurement of both the pulsar and the companion mass has been 
obtained. These values are reported in Table[2] jointly with the binary period and the 
method used to obtain the masses. 



13 



The masses of the two orbiting objects have been determined using Shapiro delay for 
six binaries out of the 11 reported in Tabled (indicated with the label "r ; s"). As we already 
discussed in §3.2^ this kind of measurement does not need any a priori hypothesis, hence 
these six objects can be used to perform six firm and independent tests on the proposed 
models. The masses for the remaining objects have been determined as follows: 



i) PSR J1012+5307 (Ivan Kerkwiik et al.lll996f ) and PSRB1957+20 (Ivan Kerkwiik et al. 



201ll ). The binary companion has been observed in phase resolved spectroscopy 



observations in the optical band. This observing mode leads to the determination 
of the companion mass function. Pulsar timing gives the pulsar mass function and 
the knowledge of these two quantities leads to the determination of the mass ratio of 
the system. The companion mass has been determined by comparing the observed 
spectroscopic properties to theoretical models for the atmosphere of HeWDs. 



201 If ). The timing 



ii) PSRs J1853+1303, J1910+1256 and J2016+1948 (IGonzalez et al. 

analysis gives a measurement of the proper motion and of the first derivative of 
the pulsar projected semi-major axis {d(ap sin i)/dt). This second quantity change s 



because of the changing orbital inclination due to proper motion (IKopeikin 



Arzoumanian et al. 



1996 . 



1996). Since this also depends on the unknown angle between the 
ascending node and proper motion directions, the most probable orbital inclination 
has been determined using Monte Carlo simulations. The companion mass is assumed 
to be that predicted by the TS99 calculation, while the pulsar mass has been obtained 
by combining the companion mass with the statistically determined orbital inclination 



and the pulsar mass function. 



In Fig. H] we plot the logarithm of the binary period versus companion mass (upper 
panels) and the pulsar mass (lower panels). We plotted all objects reported in Table|2]in the 
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left panels, while in the right ones we only show those binaries for which the masses have 
been determined from measured Shapiro delays. We marked with a triangle the position of 
PSRA in these plots to distinguish it from the field systems which are denoted by diamonds. 

In the left upper panel of Fig.0] the solid, dashed and dotted lines represent the 
predictions of the Mc — Pb correlation by TS99 (equations 20 and 21). With the exception 
of the squared point in the left bottom part of the plot, all other points well agree with 
the theoretical lines by TS99. Intriguingly the extrapolations of the results of TS99 to 
binary periods lower than 2 days are also in agreement with the observations of binaries 
with Pb < 2 days, i.e. systems whose progenitors were converging LMXBs. This means 
that the hypothesis assumed by TS99 to model the diverging LMXBs can predict rather 
well the final HeWD mass also for converging systems. The three points highlighted by a 
diamond shape represent the binaries for whic h the companion mass was obtained from 



one of the TS99 models (I Gonzalez et al 



201 ll ). so they of course fall on that line in the 



Mc — Pb plot. The squared point represents the position of PSRB1957+20, a black widow 
system. For this latter system we cannot compare the actual companion mass with the 
predictions by TS99, because of the mass loss the companion underwent after the formation 
of the MSP+WD system. Nevertheless for the binary period of this system TS99 predicts 
a companion mass ~ 0.18 M , i.e., higher than the observed one. This is consistent with 
the hypothesis that this system obeyed the TS99 Mc — Pb correlation at the epoch of the 
formation of the MSP+WD binary. 

We note that the overall agreement of the predicted Mc — Pb correlation with 
the available measureme nts has been alre ady discussed by various authors (e.g. 



van Kerkwijk et al. 



2005 and 



Lorimer 



20081 ). However these comparisons where based on a 



variety of (often) indirect methods for determining Mc, in many cases with rather large 
confidence intervals for the value for Mc- In the upper right panel of Fig.H] we display a 
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test for the Mc — Pb correlation entirely based on measures obtained via the detection 
of the Shapiro delay. The theoretical models for the Mc — Pb correlation by TS99 result 
in agreement to all 6 test binaries of our sample including PSRA, the only binary likely 
associated with a globular cluster. 

In the lower left plot of Fig.Hl eight of 11 objects show a linear correlation between 
Mp and Pb, while the remaining three, highlighted with a circle, belong to the group of 
six binaries for which the masses have been determined via the Shapiro delay detection. 
A comparison between this panel and Fig. 4-b in TS99 indicates that the theoretical 
predictions overestimate by about 0.5 M & the pulsar mass for a given orbital period in the 
entire range considered by TS99. If a real correlation exists between Mp and Pb for this 
class of binaries, it is not the one predicted by TS99 as these authors already pointed out 
in their comparison with the observed systems^. 

The solid line is a fit for all plotted points but the circled ones modelling the Mp — Pp 
correlation by using the following function: 



and we obtain A = 2.263 ± 0.005 and B = -0.466 ± 0.003 (y 2 = 0.91 with 6 degrees of 
freedom) . 

In the bottom right panel of Fig. H] the 6 binaries which represent our firm tests for 
TS99 are equally divided into two groups: three of them obey the Mp — Pb correlation 
mentioned above, while the other three do not. The binaries in this second group also 




(1) 



3 The fact that the three objects in iGonzalez et al.l (120 111 ) obey the observed Mp — Pp 
correlation can be seen as an indirect test of the results in TS99 for the Mq — Pb correlation, 
since the predictions by TS99 have been used only to obtain the companion mass. 
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have binary periods lower than the bifurcation period Pbif — 2 days theoretically found by 



Pylyser fc Savonijd (11988 



19891 ). hence their progenitors were converging LMXBs, while the 
other three have binary periods higher than P^f, so it is more likely that their progenitors 
where diverging LMXBs. 

The quantitative discrepancy between the TS99 predictions and the observed systems 
in the Mp — Pb diagram, combined with the agreement in the Mc — Pb diagram, indicates 
that in TS99 the amount of mass released by the companion is well predicted, but the 
amount of matter captured by the neutron star is overestimated. This means that the 
discrepancies between the theoretical and observed Mp — Pb diagrams result from an 
incorrect modelling by TS99 of the mechanisms (one or more) that prevent the neutron star 
from capturing all the matter released by its companion and are responsible for the mass 
loss that the binary system undergoes in its LMXB phase. Tw o mechanisms have been 



proposed, na mely the radio ejection model (IBurderi et al 



(ITauris 



20011 ) and the propeller model 



20121 ). Since our test binaries divide themselves in two separate groups, which in 



turn should have experienced different evolutionary paths for their LMXB progenitors, it is 
tempting to suggest that the mechanism responsible for the unaccreted matter is different in 
the two cases of converging and diverging LMXBs, or that in one case only one mechanism 
takes place at some point of the evolution, while in the other case both mechanisms play a 
non- negligible role. 

The position in the M P -P B plot (bottom left) of PSR J1012+5307 and PSRB1957+20 
is apparently at odds with this picture. Their orbital periods indicate that their progenitors 
where converging LMXBs, but they agree with the possible Mp — Pp correlation for binaries 
whose progenitors were diverging LMXBs. A first possibility is that the pulsar mass in 
these systems is overestimated. In fact, as reported above, Mp was obtained using indirect 
methods relying on modelling of the outer layers of the companion star. If the values 
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for M P of PSR J1012+5307 and PSRB1957+20 are confirmed by future observations, we 
could speculate that the family of converging LMXBs must be in turn subdivided in two 
groups. For the first group the mechanism(s) preventing the accretion of mass are similar 
to the diverging systems and the resulting HeWD-MSP binaries obey the same M P — P B 
correlation. The second group is instead characterised by peculiar conditions that make the 
loss of matter from the binary more efficient. Detailed numerical simulations are required 
to clarify this point. 



5. SUMMARY 

In this work we have presented a timing solution for the binary millisecond pulsar 
PSR J1910-5959A based on Parkes timing data with a data span of 12 years. Our new 
solution gives improved values for the astrometric and rotational parameters compared 
to those reported in papers I&II, and contains a new description of the orbital motion 
including the first detection of the relativistic Shapiro delay for a low-mass millisecond 
pulsar in a globular cluster. The measurement of this effect allowed us to determine reliable 
masses for the two orbiting stars. We have compared the values for the mass and the 
radius of the companion to the theoretical mass-radius relationships obtained for different 
chemical compositions, but the current uncertainties on these parameters prevent us from 
firmly establishing whether PSR J1910-5959A belongs to NGC6752 or not. Finally we 
compared the numerical results on the evolution of low-mass X-ray binary systems obtained 
by TS99, which predict the presence of a correlation between M c and Pb and between M P 
and Pb for binary systems where a millisecond pulsar orbits a low-mass helium white dwarf 
companion, to the sample of these objects for which the masses of the two stars have been 
measured using Shapiro delay. Observations confirm the M c — Pb correlation obtained by 
TS99, and show a possible M P — P B correlation for a subset of the sample which is different 
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to the one obtained by TS99. The numerical calculations by these authors overestimate the 
pulsar mass by about 0.5 M Q . These results confirm that the amount of mass lost by the 
binary during the evolution in LMXB is always larger than predicted by TS99. Moreover 
they suggest that the low-mass binary pulsars with a WD companion can be split in two 
groups, characterised by significantly different efficiencies in the mechanisms responsible for 
the loss of matter during their evolution in the LMXB phase. 

A. C, A. P., and N. D. acknowledge support for this research provided by INAF, under 
the national research grant PRIN-INAF 2010, awarded with DP 28/2011. The Parkes radio 
telescope is part of the Australia Telescope, which is funded by the Commonwealth of Australia 
for operation as a National Facility managed by CSIRO. We thank colleagues who assisted with 
the observations discussed in this paper. 
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Fig. I. — Timing residuals versus Modified Julian Date (MJD, upper panel) and orbital 
phase (lower panel) for the timing solution presented in this work. 
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Fig. 2. — Timing residuals versus orbital phase. Upper panels display unbinned residuals, 
while lower panels display residuals after being binned in 20 orbital phase bins. The left 
panels are related to the attempt of detecting the first harmonic of the Shapiro delay while 
the right panels concern the third harmonic. The solid lines in the lower panels are the 
theoretical curves for the respective harmonics. 
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Fig. 3. — Panel a) x 2 ma P i n the cos? — M c space. The cross is located at the point of 
minimum x 2 . The two solid lines correspond to the la (inner line) and 2a (outer line) regions. 
Panel b) Posterior probability density function for cost. Panel c) Posterior probability 
density function for Mq- 



Orbital period vs pulsar mass and companion mass 




M p (M ) 

Fig. 4. — Orbital period versus companion mass (upper panels) and pulsar mass (lower 
panels) plots. The left panels show the points related to all binaries in Tabled while the 
right panels show the points that represent the binary systems for which the masses have 
been determined by detecting the Shapiro delay effect only. Error bars are la uncertainties 
in the masses. If not visible, their size is smaller than the symbol. In all panels, the triangle 
represents PSRA (the only globular cluster object) diamonds represent field objects, the 
point enclosed by a square represents PSRB1957+20, the points enclosed by a diamond 
represent the three objects for which the determination of the companion mass is based on 
a model by TS99, and the circled points indicate the binaries that do not obey the linear in 
log Mp — Pb correlation (the solid line in the bottom plots). In the upper panels the solid, 
dashed and dotted lines represent the numerical results by TS99 (fig.4 in TS99 and text for 
details), while in the lower panels the solid lines represent our fit to all plotted points in the 
bottom left panel except the circled ones. 
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Table 1. Measured and derived parameters for PSR J1910-5959A. 



Parameter 


PSRJ1910-5959A 


RA ( J2000) 


19:11:42.75562(3) 


DEC (J2000) 


- 59:58:26.9029(3) 


jj, a cos S (mas yr" 1 ) 


-3.08(6) 


fig (mas yr" 1 ) 


-3.97(6) 


P(ms) 


3.26618657079054(9) 


P (io- 21 ) 


2.94703(14) 


DM (pc cm" 3 ) 


33.6998(16) 


DM1 (pc cm" 3 s" 1 ) 


0.00136(36) 


Binary model 


DD 


P orb (days) 


0.83711347691(3) 


dp sin i (lt-s) 


1.2060418(7) 


T (MJD) 


51919.206480057(55) 


M c (a) (M ) 


0.180(18) 


sin^ a ' 


>0.9994 


(deg) 


>88° 


/(Mp SR ) (M ) 


0.00268782589(5) 


(mas yr" 1 ) 


5.02(6) 


Position Angle (deg) 


217.9(7) 


Rcf. Epoch (MJD) 


51920.000 


MJD range 


51710-55911 


Number of TOAs 


1003 


Rms residuals (/is) 


5.3 



Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, 
and arcseconds. Numbers enclosed in parenthesis represent the la uncertainties on the last significant digit 
for all parameters, either measured and derived, presented in this table. The values, and relative 
uncertainties, for the parameters that describe the Shapiro delay (indicated with the apex a) have been 
obtained through a Bayesian analysis of our timing results. The proper motion's position angle is measured 

counterclockwise starting from the north direction. 
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Table 2. Masses for binary pulsars with a low mass white dwarf companion. 



PSR 



Pulsar mass (Mq) Companio mass (M©) Orbital period (days) Method 



Reference 



J0437-4715 


1.76 ±0.20 


0.254 ±0.014 


5.74 


r,s 


Verbiest et al. (2008) 


J0751+1807 


1.26 ±0.14 


0.191 ±0.015 


0.26 


r,s 


Nice et al. (2008) 


J1012+5307 


2.1 ±0.3 


0.16 ±0.02 


0.60 


optical 


van Kerkwiik et al. (1996") 


J1713+0747 


1.3 ±0.2 


0.28 ±0.03 


67.83 


r,s 


SDlaver et al. (2005) 


J1853+1303 


1.4 ±0.7 


0.33-0.37 


115.65 


TS99 


Gonzalez et al. (2011) 


B1855+09 


1.6 ±0.2 


0.270 ±0.025 


12.33 


r,s 


SDlaver (2004) 


J1909-3744 


1.438 ±0.024 


0.2038 ±0.0022 


1.53 


r,s 


Jacobv et al. (2005) 


J1910+1256 


1.6 ±0.6 


0.30-0.34 


58.47 


TS99 


Gonzalez et al. (2011) 


J1910-5959A 


1.33 ±0.11 


0.180 ±0.018 


0.84 


r,s 


this work 


B1957+20 


2.40 ±0.12 


0.035 ±0.002 


0.38 


optical 


van Kerkwiik ct al. (2011) 


J2016+1948 


1.0 ±0.5 


0.43 - 0.47 


635.04 


TS99 


Gonzalez et al. (2011) 



Uncertainties are at the la level. Pulsars are listed in right ascension order. In the column "Method" "r,s" means that the masses have 
been measured by det ecting the Shapiro delay, " optical" by analysing optical observations of the companion, TS99 by using the numerical 



results by 



Tauris fc Savoniid (|1999l ). O rbital periods have been obtained by browsing the PSRCAT pulsar catalogue 



(http://www.atnf.csiro.au/research/pulsar/psrcat, 



Manchester et al 



20051 ) . and they have been rounded to the second decimal digit since 



their uncertainties are so small that they can be considered negligible in the discussion in 
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